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Nomenclature 


B // Axial magnetic field strength-component along the line joining the contactor 

and the simulator (G) 

B j_ Transverse magnetic field strength-component perpendicular to the line joining 
the contactor and the simulator (G) 

e Magnitude of electronic charge (1.60xl0' 19 C) 

J CD Contactor discharge current (A) 

J CE Electron current emitted by contactor (A) 

J p Ion current at r e due to ions produced between r e and r B which 
flows from the potential hill region to the cathode (A) 

J SD Simulator discharge current (A) 

J SE Electron current collected by simulator (A) 

J SF Simulator filament cathode heating current (A) 

J + Ion current at r A due to ions produced between r B and r A which 
flow from the potential hill region to the downstream boundary 
(A) 

j + Current density of high energy ions flowing from the vicinity of 
the contactor to regions downstream of it (A m' 2 ) 

k Boltzmann’s constant (1.38xl0' 23 J K' 1 ) 

rh c Flowrate of neutrals supplied to contactor (seem [Xe]— standard 
cubic centimeters per minute) 

m e Mass of electron (9.11x1 0" 3 1 kg) 

m p Mass of ion (xenon: 2.18xl0' 25 kg) 

ri Neutral atom supply rate (from hollow cathode) (s' 1 ) 

n e Electron density (m' 3 ) 

n Q Neutral atom density (nrf 3 ) 

n p Density of ions on the cathode side of the potential hill [i.e. 
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those that flow toward the cathode] (m" 3 ) 

n + Density of ions on the downstream side of the potential hill [i.e. 
those that flow toward the downstream boundary] (m" 3 ) 

P 0 Ambient pressure measured far from the hollow cathode (Pa) 

R(r) Volumetric production rate of ions at radius r (s' 1 m' 3 ) 

r Radius measured from the hollow cathode (m) 

rj Radius measured from the hollow cathode (used as a dummy 

integration variable) (m) 

r A Radial position of the spherical shell at the downstream 
boundary (i.e. at the base of the potential hill) (m) 

r B Radial position of the potential hill peak (or crest) (m) 

r e Radial position of the spherical shell boundary from which 

electrons are supplied (m) 

T 0 Neutral atom temperature measured far from the hollow cathode (K) 

V Potential measured with respect io (wrt) the contactor cathode (V) 

V A Potential at r A (wrt contactor cathode) (V) 

V B Potential at r B , crest potential (wrt contactor cathode) (V) 

V c Bias supply voltage (wrt contactor cathode) (V) 

V CD Contactor discharge voltage, downstream boundary potential 
(wrt contactor cathode) (V) 

V SD Simulator discharge voltage (wrt simulator cathode) (V) 
v e Electron velocity (m s' 1 ) 

v oc Velocity of neutrals flowing from the hollow cathode (m s' 1 ) 

Z Axial position measured from the contactor cathode along the 
tank/contactor centerline (m) 

Greek symbols: 

e 0 Permittivity of free space (8.85xl0‘ 12 F m' 1 ) 

vii 


i p Solid angle of the spherical sector through which electron emission 

occurs (steradian) 

\p Q Solid angle of the spherical sector through which neutral atoms expand 
as they exit the orifice of the hollow cathode (steradian) 

<r + Electron/Neutral atom ionization cross-section (m 2 ) 

Langmuir Probe Analysis Variables and Definitions: 

A Surface area of spherical Langmuir probe- CSU probe: 3.1xl0' 5 m 2 , IFSI 
probe: 2.8xl0' 5 m 2 

E Electron energy (eV) 

E p Primary (or mono-energetic) electron energy (eV) 

F(E) Electron energy distribution function (Normalized) 

J sat Electron current flowing to a 3.1 mm dia., spherical Langmuir 
probe being held at plasma potential (A) 

n e Maxwellian electron density (cm -3 ) 

n p Primary (or mono-energetic) electron density (cm' 3 ) 

T e Maxwellian electron temperature (eV) 

V p Ep Plasma potential measured by emissive probe (V) 

V p Lp Plasma potential measured using Langmuir probe and analysis program 
developed at IFSI by Guidoni et. al. (V) 

IFSI Acronym for Istituto di Fisica dello Spazio Interplanetary of the Consiglio 
Nazionale Delle Ricerche of Italy 

CSU Acronym for Colorado State University 

LeRC Acronym for Lewis Research Center of the National Aeronautics and Space 
Administration 


ELECTRON EMISSION FROM A HOLLOW 
CATHODE-BASED PLASMA CONTACTOR 

Introduction 

Several experimental [1,2,3] and theoretical studies [4,5,6] have focused on the 
problem of controlling spacecraft electric potential with respect to an adjacent 
environment using plasma producing devices. These particular studies involved 
electrodynamic tether applications in which the plasma producing devices provide 
relatively high current capacity "contacts" to plasmas at each end of the tether. Most 
of this work concentrated on the processes that occur at the positively biased "plasma 
contactor" (i.e. the one collecting electrons from and emitting ions to the space 
plasma); and little attention was given to the negatively biased contactor that emits 
electrons. This study addresses this deficiency by focusing on the processes that occur 
near a plasma contactor emitting electrons to a simulated space plasma in a laboratory 
environment. The plasma contactor used in these experiments is a hollow cathode— a 
device derived from ion thruster neutralization applications. The hollow cathode 
plasma source is particularly well suited to charge control applications because of its 
robust construction and long lifetime characteristics, high electron emission current 
capabilities (in excess of 60 A [7]), and capacity to produce a cool, neutral plasma. 

Observations of the particles coming from a contactor emitting electrons made 
using a retarding potential analyzer (RPA) show that relatively high energy ions stream 
away along with the electrons being emitted. A mechanism is postulated in this report 
that could explain this observation. An important part of the mechanism is the high 
rate of ionization that occurs when atoms and electrons are expelled simultaneously 



through a small orifice as they are in a hollow cathode discharge [7,8]. A similar 
mechanism for the creation of high energy ions has also been proposed by investigators 
[9,10] studying various electric arcs. A plasma contactor operating in this manner 
(emitting both ions and electrons) should be well suited to spacecraft charging control 
because small changes in the potential difference between the spacecraft and the 
ambient plasma should cause the currents associated with the ion and electron flow 
between the contactor and space plasma to be altered. This would be expected in turn 
to facilitate the changes in current required to prevent the spacecraft from becoming 
either substantially positive or negative of the plasma surrounding it in a variety of 
spacecraft charging situations. This work describes recent experimental results 
obtained on a hollow cathode emitting a net electron current to a surrounding ambient 
plasma and presents a first order, one-dimensional model of the process. 

Apparatus and Procedure 

In order to study the electron emission process, the apparatus shown 
schematically in Fig. 1 was constructed. Physically this apparatus consists of two 
plasma producing devices. The one shown at the right and labeled "simulator" 
generates a simulated ionospheric plasma. The other device, shown on the left and 
labeled "plasma contactor", is the focus of this study and is biased negative relative to 
the ambient plasma to induce electron emission. Also shown are the power supplies 
and instrumentation needed to sustain and measure the characteristics of the plasmas 
produced. The simulator and contactor devices are separated by 2.7 m and are located 
within a 1.2 m dia. by 5.3 m long vacuum chamber. The contactor utilizes a 6.4 mm 
dia. hollow cathode that contains an electron emitting insert fabricated by rolling 0.013 
mm thick tantalum foil into the shape of a hollow cylinder and treating it with a low 
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a. Mechanical Schematic 


PLASMA 

CONTACTOR SIMULATOR 



Fig. 1 Experimental Apparatus 
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work function coating (containing a double carbonate [BaCC >3 , SrC0 3 ]). An orifice 
plate with a 1.0 mm dia. orifice caps the downstream end of the hollow cathode tube. 
The contactor anode is a 12 cm dia. stainless steel plate with a 1 cm O.D./5 mm I.D. 
tantalum insert near its center. The tantalum anode insert is aligned with the hollow 
cathode orifice and positioned — 2 mm downstream of it. 

Physically, the simulator resembles a ring-cusp ion source used in ion thruster 
applications [11]. Plasma is generated within it by collisions between high energy 
discharge electrons and neutral atoms. In order to increase the efficiency of this 
process, magnetic fields are used to shield anode surfaces and chamber walls against 
direct loss of discharge electrons. The ring-cusp magnetic field used in the simulator is 
induced by samarium cobalt permanent magnets. In order to ensure good coupling 
between the plasma produced within the source and the ambient plasma region, the 
device was operated without plasma extraction grids. The simulator is equipped with a 
tungsten wire cathode which is stretched diagonally across the 9 cm dia. open end of 
the source. When it is heated to thermionic temperatures it emits electrons that are 
eventually collected at the simulator body, which serves as the anode for this device as 
Fig. lb suggests. For most of the experimental results presented in this study, the 
simulator discharge current and voltage were set at 0.5 A and 40 V, respectively, and 
the simulator flowrate was set at 2.7 seem (Xe). 

Typical tests were conducted by heating the contactor hollow cathode to a 
temperature where significant thermionic electron emission could occur from the insert 
(~ 1300 K), establishing a high expellant (xenon) flowrate through it, and biasing its 
anode positive using the discharge supply to initiate a cathode-to-anode discharge. 
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Next, the desired contactor flowrate and discharge current were established; the 
contactor was biased relative to the simulator using the bias power supply; and voltage, 
current and probing instrument data were collected. The voltages and currents 
measured during typical tests are designated by the symbols shown within the circles in 
Fig. lb and defined in the nomenclature list. These quantities include the contactor 
and simulator discharge currents and voltages, the bias voltage between the contactor 
and simulator, and the contactor and simulator electron emission currents. 

The tank bias switch shown in Fig. lb was installed so the vacuum tank could 
be allowed to float relative to the contactor-simulator system or could be connected to 
the contactor cathode. Experimental results were typically not affected by the switch 
position and, consequently, it was left closed for the tests described here. The plasma 
environment produced between the contactor and the simulator was probed using the 
various instruments shown in Fig. la. They include an emissive probe, a Langmuir 
probe, and a retarding potential analyzer (RPA). The RPA consists of a cylindrical 
Faraday cage with an orifice plate at one end-the orifice hole diameter of 3 mm was 
selected to be smaller than the Debye length of the plasma in which it is typically used. 
The Faraday cage was held about 40 V below the potential of the plasma in which it 
was immersed. The probe is operated by first sighting the RPA orifice at the plasma 
contactor and then recording the ion current to the probe collector as the voltage is 
swept from 10 V below contactor cathode potential to about 100 V above it. The 
details of the current/ voltage traces obtained and their analysis and interpretation are 
discussed in Refs. [7] and [12]. 
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Experimental Observations 

Some of the phenomena observed in ground-based studies of the process of 
electron flow from a hollow cathode-based plasma contactor to a low density ambient 
plasma can be explained using the typical plasma potential profile shown in Fig. 2. In 
this case the contactor cathode, at zero potential and zero axial position, was emitting 
61 mA of electrons into an ambient background plasma located about 1 m downstream 
of the contactor. A noteworthy feature of this potential profile is the hill structure that 
develops immediately downstream of the contactor. It is postulated that this potential 
hill develops as a result of a high rate of ionization at the location of the hill. Because 
the contactor is emitting both neutrals and electrons (and both have high densities near 
the contactor), a region where the ionization rate is high can develop. Under this 
condition, electrons that cause the ionization and the electrons produced would typically 
be expected to have substantial kinetic energies after the ionization event, and they 
would be expected to leave the region quickly. However, the more massive ions would 
be left behind thereby creating a region in which the ion density is greater than the 
electron density. This net positive space-charge density region would induce a potential 
hill like the one shown in Fig. 2. It should be noted that the plasma potential data 
shown in Fig. 2 were obtained using a floating emissive probe, and these probes 
indicate potentials that fall progressively further below true plasma potential as they are 
moved into higher density plasmas [13]. Because plasma density is greatest at the 
hollow cathode orifice, the emissive probe probably indicates potentials that fall 
progressively further below true values as the cathode is approached at Z=0. Hence, 
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Fig. 2 Plasma Potential Profile on the Contactor/Vacuum Tank Centerline (Contactor 
Emitting Electrons) 


it is possible that the true crest of the hill is higher and located at a different axial 
position than the one indicated in Fig. 2. 

Downstream of the potential hill, Fig. 2 shows that the plasma potential is 
relatively constant from 15 to 60 cm. This region is called the plasma expansion 
region because the plasma density decreases in proportion to Z there and this in turn 
suggests a region of spherical plasma expansion [1]. Generally it appears that the 
plasma overexpands in this region. This is demonstrated by Langmuir probe data [1] 
which show that the plasma density at the downstream end of the plasma expansion 
region is below that of the ambient plasma region (the region of constant plasma 
potential extending beyond 100 cm in Fig. 2). It appears that the intermediate double 
layer (shown between 60 and 100 cm) enables accommodation of this difference in 
plasma densities. These phenomena and double layers in general are discussed in more 
detail by Hershkowitz [14] and the references therein. The criteria that determine the 
location, geometry and size of the intermediate double layer probably depend upon the 
ion creation and loss rates in the ambient and expanding plasma regions, the ion and 
electron current densities, and interactions with the vacuum test facility walls. 

However, the details of its characteristics have not been investigated in this study. 

A relatively low flowrate (2.3 seem [Xe]) and emission current (61 mA) were 
selected for the measurements of Fig. 2 because these operating conditions yielded a 
sufficiently low plasma density near the cathode so the potential hill could be detected 
using an emissive probe. As flowrate and/or electron emission current were increased, 
the potential hill sensed by the probe began to disappear. In order to determine if this 
was due to emissive probe inadequacy or if it indicated that the height of the potential 
hill was truly decreasing, an RPA was positioned in the plasma expansion region and 
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used to measure the energy characteristics of the ions coming from the vicinity of the 
contactor. Two typical RPA traces (recorded with the RPA positioned at Z=20 cm 
and sighted at the contactor cathode) along with their corresponding derivatives are 
shown in Fig. 3. These data were obtained with the contactor operating at a high 
flowrate (9.6 seem [Xe]) where emissive probe measurements showed no evidence of 
potential hills at either the 130 or 1000 mA electron emission levels. The RPA curve 
and corresponding derivative for the high emission current case (J CE = 1000 mA) 
indicate that two groups of ions are present. The first group induces the peak 
occurring near 15 V in the lower plot and represents low energy, thermal ions present 
in the expanding plasma. The second group which exhibits a greater energy spread is 
present as the tail on the solid curve extending from 20 to 100 V (Fig. 3b). It is 
postulated that the high energy ions associated with this tail are created on a potential 
hill located near the contactor cathode and that they flow from there to the RPA where 
they are detected. The RPA data corresponding to contactor operation at a low 
electron emission current of 130 mA displays only one low energy group of ions. This 
suggests that no potential hill forms at this operating condition. 

There are other differences between the plasmas measured in the expanding 
plasma region at 130 and 1000 mA emission currents and one of these, the difference 
in electron energy distribution functions sensed by a Langmuir probe, is illustrated in 
Fig. 4. At a high emission current, the solid curve suggests that two electron groups 
exist. One group, associated with the lower energy peak, probably represents the 
thermal electrons present in the expanding plasma. The other, higher energy group is 
associated with electrons that have been accelerated from the contactor cathode through 
the potential hill region and into the expanding plasma without experiencing many 
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RPA COLLECTOR POTENTIAL (V) 


a. Current/Voltage Traces 



RPA COLLECTOR POTENTIAL (V) 
b. Derivatives of Current/ Voltage Traces 

Fig. 3 Retarding Potential Analyzer Data Measured in the Expanding Plasma Region (High 
Contactor Flowrate Condition) 
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Fig. 4 Electron Energy Distributions Measured in the Expanding Plasma Region (High 
Contactor Flowrate Condition) 



energy dissipating collisions. However, the electron distribution function 
corresponding to the low emission current condition (J CE =130 mA) indicates only one, 
low energy group of electrons is present. Thus, Figs. 3 and 4 show that both ions and 
electrons in the expanding plasma region exhibit distribution functions that have 
thermal and high energy components at a high emission current, while only the thermal 
component is present at a lower electron emission current. Note that the electron 
energy distribution functions shown in Fig. 4 have been normalized with respect to 
their maxima. They were computed from the second derivatives of spherical Langmuir 
probe traces recorded digitally using a Keithley 617 programmable electrometer and 
plasma potential data measured using the emissive probe. The derivatives were 
obtained by finding the discrete Fourier series representation of the Langmuir probe 
trace, solving for the analytical derivative of this series and applying a convergence 
filter [15]. The convergence filter artificially smooths the data and tends to spread out 
quickly varying features of the Langmuir probe trace. Smoothing errors introduced by 
the data reduction procedure and plasma potential measurement errors yield distribution 
functions (Fig. 4) that reflect only the qualitative differences between the plasmas 
observed at the 130 and 1000 mA operating conditions. Appendix A contains a more 
detailed explanation of the procedures used to compute electron energy distribution 
functions, and Appendix B contains an explanation of the technique used to obtain 
plasma potential using a floating emissive probe. 

It is also of interest to examine the effect of electron emission current on the 
axial profiles of the high energy ion current density and the random electron current. 
These profiles have been measured over a range of electron emission currents and the 
results are shown in Fig. 5. The data in this figure correspond to a lower flowrate 
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ELECTRON SATURATION HIGH ENERGY ION CURRENT 

CURRENT [J sat ] (yUA) DENSITY [jj (yUA, cm" 2 ) 


^CD — ^ A 



a. Effect of Electron Emission Current on High Energy Ion Current Density Profiles 



b. Effect of Electron Emission Current on Langmuir Probe Electron Saturation 
Current Profiles 


Fig. 5 Experimental Observations Suggesting Spherical Expansion in the Region Downstream 
of the Potential Hill Structure . 0 



(4.1 seem [Xe]) than those of Figs. 3 and 4. At this lower flowrate, high energy ions 
were detected flowing from the contactor at all four of the electron emission current 
levels shown (i.e. at J CE = 126, 500, 1000 and 1500 mA). The lines drawn on the two 
plots shown in Fig. 5 correspond to an inverse square dependence with axial position. 

The high energy ion current density is shown to follow the inverse square 
dependence (Fig. 5a) and this suggests that the high energy ions are indeed expanding 
spherically from their point of creation. On the other hand, the electron saturation 
current in the expanding plasma region decreases less rapidly than the inverse square of 
distance— especially at large values of axial position as shown in Fig. 5b. The electron 
saturation current [J sat ] is defined to be the electron current flowing to the Langmuir 
probe held at plasma potential. Hence this current includes both high energy (beam) 
electron and thermal electron components. It is considered likely that the data of Fig. 
5b do not always follow the Z‘ 2 dependence either because the beam electrons are 
being thermalized or their presence is being masked by thermal electrons in the plasma 
expansion region. 

Both the high energy ion and directed electron populations present in the 
expanding plasma can contribute to instabilities and cause the expanding plasma to be 
noisy or turbulent. A coarse measure of the turbulent intensity (the fraction of the 
energy in the expanding plasma that is in the form of turbulent electrostatic 
fluctuations) is the square of the ratio of the rms density fluctuations to the mean 
plasma density. This density ratio can be measured qualitatively in the low density 
expanding plasma by monitoring the current to a Langmuir probe when it is held near 
plasma potential and recording the rms noise amplitude/mean current ratio. Figure 6 
shows rms-to-mean current ratio versus axial position data corresponding to the data 
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Fig. 6 Langmuir Probe Noise Levels at High and Low Electron Emission Currents (High 
Contactor Flowrate Condition) 


presented in Figs. 3 and 4. The data for the 1000 mA operating condition suggest that 
the plasma is very noisy near the plasma contactor (turbulent intensity ~ [0.32] 2 = 
10%) and less noisy (~ 2%) at larger axial positions, and the opposite trend is 
indicated for the 130 mA operating condition. The noise levels at 33 cm are shown to 
be comparable thereby suggesting that phenomena occurring in the ambient plasma 
region determine the noise level at axial positions greater than 30 cm. 

Theoretical Development 

In order to gain some understanding of the potential hills that have been 
measured at low emission currents and postulated at higher ones, a simple model of the 
electron emission process has been developed. Figure 7 is a sketch of the spherical 
geometry associated with a hypothetical potential profile that shows the electron 
emission current flowing from the hollow cathode through the potential hill to a 
downstream boundary. As these electrons leave the electron source surface, they are 
accelerated up the potential hill and they gain sufficient energy to ionize neutral atoms. 
The resulting ions will flow down the hill from the point where they were produced. 
Ions produced to the left of the crest potential shown in Fig. 7 will flow to the cathode 
and those produced to the right of it will flow to the downstream boundary. The 
streaming electrons, which accelerate to the crest and then decelerate after they pass it, 
will still have substantial kinetic energies as they pass the downstream boundary. 

The approach used here will be to write equations that describe the electron and 
ion densities throughout the region between the electron source and the downstream 
boundary and then apply Poisson’s equation to solve for the associated potential profile. 
Because the electron and ion densities depend upon the potential profile, however, an 
iterative solution technique must be applied to accomplish this and obtain the 
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Fig. 7 Spherical Model Diagram of the Electron Emission Process 
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steady-state, self-consistent solution for the density and potential profiles. This model 
of the electron emission process will be presented in terms of two sets of equations. 
One set will pertain to radial locations between the cathode and the potential peak (i.e. 
on the cathode side of the potential hill). The other set will pertain to radial locations 
between the potential hill and the downstream boundary (i.e. on the downstream 
boundary side). 

The Cathode Side (r e < r < r B ) 

The electron density at any point in this region can be described approximately 
by assuming conservation of electron energy and current, 
i.e. 


I m e [v e (r)] 2 = e V(r) 


( 1 ) 


and 

J CE = e n e( r ) f r v e( r ) 


■ ( 2 ) 


Combining these equations and solving for the electron density gives 


n e (r) 


J CE m e 


• ( 3 ) 


This expression for the electron density is only approximately correct because it ignores 
both electrons which are produced in ionization events and the effects of energy 
removal from the electron group due to ionization and other inelastic collisions. 
Neglecting these effects to make the problem more tractable limits its direct 
applicability to the case where the inelastic collision rate expressed as a current is small 
compared to the electron emission current. It is assumed that some mechanism for 
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removing low energy electrons produced via ionization from the potential hill region is 
active. Although this mechanism is not defined, it is noted that the current of these 
electrons should be small so a negligible fraction of the kinetic power in the streaming 
electrons would be required to remove them through elastic collisions. 

The rate of ion generation per unit volume [R(r)] at radius r is given by 

R(r) = n e {r) n Q (r ) <r + (v e ) v e (r) • (4) 

The ionization cross-sections of Rapp and Englander-Golden [16] were used to calculate 
ion generation rates. In addition, the neutral atom density n Q appearing in this equation 
is assumed to be the sum of the neutral densities due to the neutral flow from the 
hollow cathode (assumed to expand spherically from the orifice) and the background 
neutral density in the vacuum environment of the test. Specifically, the density n 0 at 
radius r was approximated by 


n 0 (r) 



• (5) 


The density of ions at a radius r is determined by summing the contributions of all ions 
produced at radii of greater potential. Each of these ions will be accelerated from their 
point of creation rj to the radius of interest r. Hence, the contribution to the density of 
ions at a radial location r (for the region r e < r < r B ) due to ions generated with a 
negligible initial velocity in a differential volume near r t is 


dn p (r) = r~ 2 


R( r i) r\ dr j 

Wi) 


• ( 6 ) 
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and the velocity of the ions created at r 1 once they reach r is given by 


vOtf 


2e[V(r l )-V(r)] 


• (7) 


The overall ion density at any radius r on the cathode side of the hill is now found by 
integrating the differential density dn p from r to r B . This yields 



*0i) r \ d r \ 

Wo 


• ( 8 ) 


Combining Eqs. 3 through 8 and simplifying gives 


n D {r) = r 


-2 J CE 


■ji * 


,3/2 


*A v oc 


k T n 


[V(ri)-V(r)] 


- 1/2 


dr x 


(9) 


The electron and ion densities determined using Eqs. 3 and 9 can now be combined 
with Poisson’s equation to describe the variation in plasma potential on the cathode side 
of the potential hill. Assuming spherical symmetry, Poisson’s equation is 


d 2 V 


li Z-~ w 

r dr ' 



• ( 10 ) 


In order to apply the equations just developed, it is necessary to develop the 
equations describing conditions on the opposite (downstream boundary) side of the 
potential hill so all of them can be solved simultaneously. 
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The Downstream Boundary Side (r B < r < r A ) 

Under the assumptions of this development, the equation that describes the 
electron density in the region between the potential crest and the downstream boundary 
is the same as the one developed for application upstream of the potential crest, namely 
Eq. 3. The ion density expression is obtained by repeating the logical sequence used to 
derive Eq. 9. It is found to differ from Eq. 9 only in the order of the integration, 
hence 

[Wi )-W)Y m ( 11 ) 

Note that Eq. 11 shows an inverse square dependence with position and a linear 
dependence with emission current. This is in qualitative agreement with the functional 
dependencies indicated by the experimental data of Fig. 5a. 

Equations 3, 9, 10, and 1 1 represent a relatively simple model of the electron 
emission process. They were solved by first dividing up the region r e to r A using 
small, evenly-spaced node points. Next, the derivatives in Eq. 10 were approximated 
using finite-difference expressions. This allowed algebraic equations arranged in 
matrix form to be written for the potential at each node point. Electron and ion 
densities were then calculated at each node point using Eqs. 3, 9 and 11 and an initial 
estimate of the potential variation through the potential hill region. The procedure of 
solving for the densities and then the potentials at each node was repeated many times 
until the potential profile stabilized. 

It should be noted that the solution procedure treats 1) the electron source 
location r e , 2) the downstream boundary location r A , 3) the solid angles \p and \ p Q , and 


n + (r) = r 


-2 J CE 


ft * 


,3/2 


r B 


*o r \ v oc 


k T„ 
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4) the potential at the downstream boundary V A as parameters. The electron source 
and downstream boundary locations are, however, not free parameters. The values of 
r e and r A are established physically by the requirement that the electric fields be zero at 
these locations (i.e. the space-charge limited condition applies). It was postulated that 
the other parameters, namely the downstream boundary potential V A , and the solid 
angles $ and \ p Q were influenced by such factors as the cathode orifice size, the anode 
configuration, and the plasma conditions beyond the downstream boundary and they 
were treated as free parameters. It is believed that an energy balance analysis could be 
used to find the downstream boundary potential, but this was not accomplished in this 
preliminary study. In order to apply the model and compare its predictions to 
experimental observations, V A was set at the experimentally measured potential in the 
expanding plasma region (typically measured at a radius of 20 cm) for each electron 
emission operating condition studied. The solid angles and \p 0 were arbitrarily set to 
2 t (i.e. hemispherical geometry). There are other parameters appearing in equations 
3, 9, 10, and 11 that are not determined explicitly through the analysis (e.g. P 0 and 
T 0 ), but they were controlled during the experiment and unique values could be 
assigned to them. 

Numerical Example 

When Eqs. 3, 9, 10, and 11 were solved for the case of an emission current of 
1 A and values of the parameters given in Table 1, the theoretical potential profile 
shown in Fig. 8a was computed. By forcing the boundary electric fields to be zero, 
the electron source and downstream boundaries were found to be located at 4.6 and 
14.4 mm, respectively, and the crest potential of 153 V was located at 7.4 mm. This 
large potential was caused by a net positive (ion) charge density in the region between 
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Table 1 Numerical Example Data Set 


Inputs 

Outputs 

^CE = 1*0 A 

Fig. 8 

m c = 4.1 seem (Xe) 
(or n = 1.72xlO l8 s' 1 ) 

J + = 1.94 mA 


J = 2.16 mA 

P 0 = 5.0X10- 6 Ton 
(or 6.7X10" 4 Pa) 

V 

r e = 4.6 mm 

T 0 = 300 K 

r B = 7.4 mm 

v oc = 458 m s" 1 

r A = 14.4 mm 

II 

to 

V B = 153 V 

tN 

II 

o 


V A = 22 V 
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5 and 11.5 mm as shown in Fig. 8b. The neutral atom density variation throughout the 
potential hill region is shown in Fig. 8c. When this information was combined with 
the data shown in Figs. 8a and 8b the ion production rate per unit volume was 
calculated and it is plotted in Fig. 8d. By integrating the volumetric ion production 
rate over the entire potential hill region volume, the total ion current flowing from this 
region was calculated to be 4.10 mA. An ion current of about 1.94 mA (J + ) was 
found to be flowing from the potential hill region to the downstream boundary and 2.16 
mA (J p ) was found to be flowing to the electron source boundary. These ion creation 
rates (expressed as currents) are small compared to the emission current. This suggests 
that the assumptions made in deriving this model are probably valid and that very little 
power should be required from the streaming electrons to remove low energy electrons 
(resulting from inelastic collisions) from the potential hill as quickly as they are 
produced. 

Comparison of Theory and Experiment 

The procedures used to obtain the numerical results given in Table 1 and Fig. 8 
were applied to obtain additional solutions over ranges of electron emission currents 
and flowrates. The effect of electron emission current and flowrate on the current 
density and maximum energy of ions flowing away from the hollow cathode discharge 
were also measured using the RPA described previously. The experimentally measured 
and theoretically predicted effect of emission current on these quantities are compared 
in Fig. 9 under conditions where the RPA was positioned 20 cm downstream of the 
contactor. Figure 9a shows the high energy ion current density increasing with 
electron emission current, at a smaller slope than the "theoretical" curve. The 
theoretical predictions of high energy ion current density were made by first finding the 
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ion current emitted from the potential hill region to the downstream boundary for the 
particular electron emission current as explained in the numerical example of Table 1 
and Fig. 8. Next, this current was divided by the surface area at a radius of 20 cm 
(i.e. t^r 2 * 2ir[20] 2 cm 2 « 2500 cm 2 ) to obtain the current density at this location. 
Although the two curves shown in Fig 9a do not coincide, the agreement between the 
experiment and numerical model is considered to be good considering the assumptions 
made in the model. Uncertainties in experimental conditions as well as in the 
ionization cross sections could easily cause the level of error indicated in Fig. 9a. It is 
noted that better agreement could be artificially generated in this simple one- 
dimensional model by adjusting the solid angle yp with each electron emission current. 
However, it is felt that two-dimensional (or possibly three-dimensional) effects 
probably determine the subtle trends in the experimental data so attempts to adjust \p to 
obtain better agreement would not be justified. 

Figure 9b contains a comparison of experimentally and theoretically 
determined crest potentials. Again, relatively good agreement and a similar trend for 
the crest potential to increase with electron emission current for both curves is shown. 
The actual positions of r e , r B , and r A at the electron emission levels corresponding to 
Fig. 9 are shown in Fig. 10. The most notable trend in this figure is that larger values 
of r A correspond to smaller values of electron emission current. Together with Fig. 

9b, this suggests that not only are crest potentials greater at higher electron emission 
currents, but electric field strengths are also higher. 

The effect of contactor flowrate on the experimentally measured and 
theoretically predicted high energy ion current and crest potential are shown in Fig. 11. 
The theoretical predictions (triangular and solid circular data points) and experimental 
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Fig. 10 Computed Values of Downstream, Crest, and Electron Source Boundaries as 
Functions of Electron Emission Current 


28 



0.8 


^CD =b O' 1 A 
V CD = 8 TO 14 V_ 

P Q = 3 TO 7 xIO 6 Torr 
Z=* 20 cm 



a. Ion Current Density 
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Fig. 11 Computed Effects of Contactor Flowrate 
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measurements (circular and square data points) of high energy ion current density 
corresponding to the 130 and 1000 mA electron emission levels are shown to exhibit 
comparable magnitudes in Fig. 11a. The theoretical crest potential data for the 1000 
mA condition shown in Fig. lib also show good agreement with experiments, and the 
predicted and measured crest potentials show a similar trend (both decrease with 
flowrate). However, the predicted and measured crest potentials corresponding to the 
130 mA condition do not show the same trends. Note that current density 
measurements made at an electron emission current of 130 mA (shown in Fig. 11a) 
indicate, at a contactor flowrate of 9.6 seem, that no high energy ions are present. At 
this high flowrate and low electron emission level, apparently no potential hill structure 
is needed to assist electron emission from the high density plasma near the contactor 
hollow cathode. Numerical modelling of the 130 mA condition was difficult to 
perform at the higher flowrates and, in order to obtain a steady solution, the 
downstream potential V A had to be artificially increased to — 14 V (from actual 
measured values of 8 to 12 V) to realize a convergent solution. This action probably 
caused the relatively constant value of crest potential (about 55 V as shown in Fig. 

1 lb) to be exhibited at high flowrates. For the convenience of the reader, the predicted 
values of r e , r B , and r A corresponding to Fig. 11 are shown in Fig. 12. It shows that 
higher flowrates induce larger radii and these values appear to grow linearly with 
flowrate. 

The level of agreement between the model and experiment shown in Figs. 9 
and 11 is considered to be good. In addition, the rate of ion production (expressed as a 
current) in the potential hill region is small compared to the electron emission current 
for all of the comparisons made in these figures. Hence, the rate of low energy 
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Fig. 12 Computed Values of Downstream, Crest and Electron Source Boundaries as 
Functions of Contactor Flowrate 
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(secondary) electron production is also small and it should be possible to remove them 
as they are produced using a small fraction of the power in the streaming electron 
beam. Recall their removal via an undefined mechanism is assumed in the model. 

Conclusions 

Experimental observations of a hollow cathode-based plasma contactor emitting 
electrons to a ambient plasma suggest that a potential hill structure develops close to 
the contactor cathode. It is postulated that the potential hill is created by a region of 
positive space charge, and ions produced in this region can gain substantial energies as 
they are accelerated away from the region. By measuring the energies of these ions, 
the height of the potential hill can be inferred. In general, an increase in contactor 
flowrate tends to reduce the potential at the crest of the hill, while larger emission 
current levels tend to increase it. 

A simple model that reflects the effects of ionization, ion and electron 
acceleration and the space-charge induced by the ions and electrons describes the 
essential features derived from experimental observations of hollow cathodes emitting 
electrons. Specifically, it yields magnitudes of potential hill height and current density 
of ions flowing from the potential hill that agree with experimental results. Further, 
the predicted effects of electron emission current and contactor flowrate on these 
features agree with experimentally observed trends. It is noted that the total current of 
ions emitted to the expanding plasma is estimated to be small compared to the electron 
emission current (i.e. typically less than 0.2 %). This suggests that only a small 
fraction of the electrons flowing from the contactor to the expanding plasma interact 
with and possibly ionize neutral atoms while they stream through the potential hill 
region. 
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CSU AND CNR-IFSI BI-LATERAL 
EXPERIMENTS ON PLASMA CONTACTORS 

Introduction 

Experiments have been conducted independently at Colorado State University 
[1], NASA Lewis Research Center [3], and the Italian Institute of Interplanetary Space 
Physics [2,19] on the electron collection characteristics of hollow cathode-based plasma 
contactors. The most basic experiments involve biasing a plasma contactor and its 
associated plasma cloud (or high density plasma plume) positive with respect to an 
ambient plasma and measuring the current which is collected under this applied 
voltage. Typically, most of the voltage drop between the contactor (and its plasma 
cloud) and the ambient plasma develops across a double layer. Previous experiments 
conducted at CSU, NASA Lewis, and IFSI have detected double layers, but some 
differences have been noted. In order to address these differences, a bi-lateral research 
program was initiated between CSU and IFSI. This report, which describes work 
conducted under this program, presents the results of recent experiments conducted in 
the IFSI facility on the electron collection mode of the plasma contacting process. The 
goals of the IFSI experiments were to 1) investigate the effects of magnetic field on the 
electron collection process, 2) compare and calibrate Langmuir and emissive probe 
plasma diagnostic techniques, and 3) operate a hot filament cathode-based ion source as 
an ambient plasma generator. 

Apparatus and Procedures 

In order to study the electron collection process, the apparatus shown 
schematically in Figs. 13 and 14 was set up at IFSI. Physically this apparatus is 


33 



34 


Fig. 13 IFSI Vacuum and Probing System Schematic 
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Fig. 14 IFSI Electrical Schematic 




similar to the one described previously in this report. The simulator and contactor 
devices are separated by 2.7 m and are located within a 2 m dia. by 4 m long, stainless 
steel vacuum chamber. The contactor utilizes a 6.4 mm dia, orificed hollow cathode 
which houses an electron emitting Ta foil insert. The hollow cathode orifice is 0.38 
mm in diameter and the contactor anode is a 12 cm dia. stainless steel plate with a 1 
cm dia. tantalum insert that has a 5 mm dia. orifice in it. The anode plate and its 
tantalum insert are located concentric with the hollow cathode centerline at a plane —2 
mm downstream of the cathode orifice plate. 

Typical tests were conducted by heating the contactor hollow cathode to a 
temperature where significant thermionic electron emission could occur from its insert 
(— 1300 K), establishing an expellant (xenon) flowrate through it, and biasing its anode 
positive using the discharge supply to initiate a cathode-to-anode discharge. The 
simulator, a hot-filament equipped ion source shown in Fig. 15 and described in Ref. 

12, was also started. Next, the desired contactor and simulator flowrates and discharge 
current levels were established; the contactor was biased relative to the simulator using 
the bias power supply; and voltage, current and probing instrument data were collected. 
The tests described in this attachment were performed at a contactor flowrate and 
discharge current of 1.2 seem (Xe) and 0.6 A, respectively. This operating condition 
induced a 30 V discharge voltage, and a background pressure of 16xl0' 6 Torr. The 
simulator flowrate was about 2 seem (Xe), and its discharge current and voltage were 
set to 0.6 A and 54 V, respectively. 

The plasma environment produced between the contactor and the simulator was 
probed using the instruments shown in Fig. 13. The emissive probes were used to 
measure plasma potential and they were compared to Langmuir probe plasma potential 
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Fig. 15 Hot Filament Cathode-based Simulator 
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measurements. The nine Langmuir and two emissive probes were fixed to a movable 
platform that could position them on a line that intersected the centerline and was 
parallel to the plane of the contactor anode at axial positions varying from 1 to 265 cm 
(measured from the contactor anode). Figure 16 is a perspective view of the probes 
showing their locations in relation to the plasma contactor. The center Langmuir probe 
was positioned directly in front of the plasma contactor orifice, and the two emissive 
probes were positioned 3 cm from the axis on each side. The nine Langmuir probes 
were equally spaced at intervals of 12 cm. 

The magnetic field present in the region between the contactor and simulator 
could be controlled in both magnitude and direction by large Helmholtz coils which 
encircled the stainless steel vacuum tank. In order to study its effect on the electron 
collection process, various magnetic field configurations were imposed. They included 
1) a zero magnetic field, the geo-magnetic field was nulled; 2) axial fields, those 
directed along the axis joining the contactor and simulator; and 3) transverse fields, 
those aligned perpendicular to the axis joining the contactor and simulator. 

Results and Discussion 

Figures 17 and 18 show typical plasma property data that were either 
measured directly or computed from measurements made in the IFSI laboratory. 

Figure 17 corresponds to a condition where the Earth’s magnetic field was nulled and 
100 mA of electrons were being collected by the contactor from the ambient plasma. 
The emissive probe data in the upper plot of Fig. 17 show the potential variation 
through the double layer that develops between the contactor plasma cloud and the 
ambient plasma. It reveals a double layer across which a 27 V potential difference 
develops in the region between Z = 15 and 30 cm. Two energy distribution function 
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Fig. 16 Emissive and Langmuir Probe Arrangement 





plots derived from Langmuir probe data that describe the electrons at axial positions on 
either side of the double layer (10 and 49 cm) are shown adjacent to the plasma 
potential profile. The downstream distribution function (at 49 cm) suggests that only 
one group of electrons is present in the ambient plasma and they appear to be 
Maxwellian. On the other hand, the upstream distribution function (at 10 cm) shows a 
high energy group in addition to the Maxwellian one. It is believed that electrons 
accelerated from the ambient plasma through the double layer and into the upstream 
contactor plasma cloud region, are the ones that induce the non-thermal (high energy) 
hump centered near 30 eV. 

The middle and bottom plots in Fig. 17 contain electron density and electron 
temperature data obtained from Langmuir probe measurements made in the contactor 
plasma cloud and ambient plasma regions. The Maxwellian electron group density is 
shown to be quite high near the contactor and to decrease rapidly with increasing axial 
position. The primary electron density is lower and it remains relatively constant 
within the contactor plasma cloud. It is noted that the densities of the Maxwellian and 
primary electrons in the contactor plasma cloud were calculated using the CSU 
Langmuir probe analysis program, which is based on the assumption that two electron 
groups are present in the plasma. One group is modelled as Maxwellian, while the 
other one (the primary group) is assumed isotropic and mono-energetic. The program 
solves for the Maxwellian group temperature and density and the primary group energy 
and density by using a non-linear, least-squares curve-fit to the portion of the Langmuir 
probe data in the electron retarding region. The properties of the electrons were also 
analyzed using a computer program developed by U. Guidoni, et. al. at IFSI. 

Although two separate Maxwellian groups can be used in the IFSI program to model 
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the plasma, the program was used to analyze only the colder Maxwellian electron 
groups in the contactor plasma cloud and ambient plasmas. Results obtained using 
these two analysis techniques will be compared shortly. 

The data in the middle plot of Fig. 17 show the Maxwellian electron density in 
the ambient plasma region remains relatively constant with axial position. 

Temperatures computed from Langmuir probe measurements for the Maxwellian 
electron groups in the ambient plasma and contactor cloud regions are shown in the 
lowest plot. It reveals a very low electron temperature in the ambient plasma (about 
0.7 eV) and an electron temperature that rises in the contactor cloud from 1.8 eV close 
to the contactor to 5 eV at a position near the double layer. This rise in electron 
temperature is consistent with observations made in other double layer tests [14]. It 
suggests electron heating is occurring near the double-layer boundary possibly as a 
result of turbulent interactions between the high energy and Maxwellian electrons 
present in the contactor plasma cloud. Indications of higher electron temperatures 
could also be due to an error in the Langmuir probe analysis program that becomes 
significant when the high energy electron signal begins to dominate the colder electron 
group signal at locations close to the double layer boundary. 

Figure 18 contains data which correspond to a transverse magnetic field with a 
magnitude of 0.35 G. It is organized similar to Fig. 17 and it shows many of the 
trends attributed to Fig. 17. However, one important difference between Figs. 17 and 
18 is the position of the double layer. It is located further downstream of the contactor 
(between 26 and 35 cm) when the transverse field is applied. Presently, the processes 
which cause the contactor cloud to enlarge in the axial direction when a transverse field 
is applied are not understood. 
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Extensive Langmuir probe data were collected at IFSI (about 1500 traces!) for 
several electron collection currents and magnetic field configurations. A comparison of 
the results obtained using both the IFSI and CSU analysis programs for Langmuir 
probe traces taken on the centerline is included in Tables 2 and 3. These tables also 
include emissive probe measurements of plasma potential (V p EP ), which are typically 
used in the CSU analysis routine but not in the IFSI one. These plasma potential data 
are seen to agree with values obtained from Langmuir probe measurements found using 
the IFSI program (V p LP ). It is noted, however, that a significant difference tends to 
develop between plasma potentials determined from the Langmuir and emissive probe 
measurements at small values of electron collection current— see the top data set of 
Table 2, for example. This error is believed to be caused by the fact that the emissive 
and Langmuir probes do not measure potential at the same location (they are separated 
by 3 cm). When the electron collection current is large, the plasma contactor cloud is 
also large and both probes tend to lie within it when they are close to the contactor. At 
low emission currents, however, the emissive probe may lie outside of what tends to be 
a small contactor plasma cloud. It is also noted that it is probably impossible to 
interpret Langmuir probe data taken within and very near the double-layer region. It is 
because of this that the data of Figs. 17 and 18 and Tables 2 and 3 do not show 
electron densities and temperatures/energies in this region. 

In general, the agreement between the IFSI and CSU Langmuir probe data 
analysis programs is considered to be good. The largest differences are observed 
between Maxwellian temperature results in the contactor plasma cloud region. A 
contributor to this difference was the large voltage step (IV) used in the acquisition of 
the Langmuir probe data. As mentioned previously, another difficulty encountered 
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Table 2 Centerline Plasma Property Comparisons for Nulled Magnetic Field 
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40 

6.0 

6.2 

1.1x10® 

0.95 

— 

— 

5.6xl0 7 

1.2 

vd043 

43 

5.8 

6.2 

1.0x10® 

0.95 

— 

— 

6.3xl0 7 

l.l 

vd046 

46 

5.8 

6.3 

7.9xl0 7 

1.1 

— 

— 

6.4xl0 7 

1.2 

vd049 

49 

5.8 

6.3 

7.6x1 0 7 

1.1 

— 

— 

6.5xl0 7 

1.2 

vd055 

55 

5.8 

6.3 

■SSEKI 

1.0 

— 

— 

6.8xl0 7 

1.1 

vd061 

61 

5.8 

6.3 


1.0 

— 

— 

6.9xl0 7 

l.l 


45 

































































































































































































































Table 3 Centerline Plasma Property Comparisons (Zero, Transverse, and Axial Magnetic Fields) 


Jqj? * -100 mA 
B x »0G 
B, - 0 G 

] < CSU ANALYSIS > | <- IFSI ANALYSIS -> j 


Ref.# 

Z 

V P,EP 

V P,LP 

n e 

T, 

"p 

H 


T, 


(cm) 

(V) 

(V) 

(cm* 3 ) 

(eV) 

(cm* 3 ) 

(cV) 


(eV) 

re 001 

1 

31 

35 

6.5xl0 9 

1.8 

3.2xl0 7 

33 

6.5xl0 9 

1.9 

rc004 

4 

31 

33 

1.6xl0 9 

1.7 

2.2xl0 7 

32 

1.3xl0 9 

2.2 

rc007 

7 

31 

34 

6.5x10 s 

2.1 

2.8xl0 7 

32 

4.3xl0 8 

3.7 

vc013 

13 

30 

36 

2.0xl0 8 

4.3 

3.9xl0 7 

35 

1.9xl0 8 

8.7 

vc016 

16 

26 

34 

1.4xl0 8 

4.9 

3.4xl0 7 

33 

1.5xl0 8 

9.2 

vc034 

34 

5.2 

5.5 

7.2xl0 7 

0.72 

— 

— 

7.6xl0 7 

0.86 

vc037 

37 

5.1 

5.5 

6.7xl0 7 

0.72 

— 

— 

9.5xl0 7 

0.96 

vc 049 

49 

5.1 

5.5 

1.4xl0 8 

0.64 

— 

— 

7.4xl0 7 

0.91 

vc061 

61 

5.1 

5.5 

1.8xl0 8 

0.68 

— 


8.6xl0 7 

0.89 


J r»p = -100 mA 
B x - 0.35 G 
B, - 0G 

J< CSU ANALYSIS > j <- IFSI ANALYSIS -> j 


Ref. # 

Z 

V P,EP 

V P,LP 


T, 

d p 

wm 

n e 

T e 

re 001 

1 

33 

36 

8.4xl0 9 

1.5 

1.6xl0 7 

36 

6.0X10 9 

1.7 

re 004 

4 

34 

36 

2.6xl0 9 

1.4 

1.6xl0 7 

35 

1.2XI0 9 

2.1 

ve007 

7 

34 

37 

6.5xl0 8 

2.0 

1.8xl0 7 

35 

4.8x10 s 

3.0 

ve013 

13 

33 

37 

2.4xl0 8 

2.7 

I.4xl0 7 

35 

1.7x10 s 

4.8 

ve019 

19 

33 

37 

1.6xl0 8 

3.1 

1.5xl0 7 

35 

1.2x10? 

6.1 

ve025 

25 

32 

37 

8.6xl0 7 

3.8 

1.8XI0 7 

36 

7.0xl0 7 

9.6 

ve046 

46 

6.0 

5.5 

4.9xI0 7 

0.99 

— 

— 

2.0xl0 7 

1.3 

ve052 

52 

4.8 

5.1 

2.9xl0 7 

1.0 

— 

— 

2.1xl0 7 

l.l 

ve058 

58 

4.7 

5.6 

6.7xl0 8 

0.95 

— 

— 

3.6xl0 7 

1.1 

ve061 

61 

4.6 

5.5 

1.2xl0 8 

0.82 

— 

— 

5.0xl0 7 

1.1 


= -100 mA 
B ± - 0 G 
B, = 0.35 G 

j < CSU ANALYSIS > | < - IFSI ANALYSIS -> j 


Ref. # 

Z 

V P,EP 

V P,LP 


T, 

“p 

WM 

n e 

T, 

rkOOl 

1 

33 

33 

5.2xlC? 

1.4 

2.1xl0 7 

31 

3.2X10 9 

1.8 

vk004 

4 

34 

33 

8.8x10? 

1.8 

1.3xl0 7 

31 

6.0x10 s 

2.5 

vk007 

7 

34 

33 

3.4x10? 

1.9 

l.7xl0 7 

32 

2.5x10 s 

3.5 

vk013 

13 

34 

34 

3.4x10? 

2.4 

2.9xl0 7 

33 

1.6x10 s 

5.9 

vk016 

16 

30 

34 

2.7x10 s 

2.5 

4.1xl0 7 

34 

1.3x10 s 

9.5 

vk034 

34 

5.0 

6.0 

4.6xl0 7 

1.7 

— 

— 

5.7xl0 7 

1.2 

vk037 

37 

4.9 

6.0 

8.8xl0 7 

1.5 

— 

— 

6.3xl0 7 

1.2 

vk040 

40 

4.8 

5.9 

1.1x10? 

1.3 

— 

— 

6.5xl0 7 

1.2 

vk043 

43 

4.8 

5.9 

1.1x10? 

1.3 

— 

— 

6.7xl0 7 

1.2 
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when attempting to evaluate the properties of the colder electron group present in the 
contactor plasma cloud was the perturbation and masking of this signal by the high 
energy electron group signal. 

The data of Table 2 were measured when the magnetic field strength in the 
region between the contactor and simulator was reduced to zero. The plasma property 
measurements that yielded these data were made along the tank centerline at electron 
collection currents of 50, 100, and 150 mA. Table 3 contains data obtained at an 
electron collection current of 100 mA in magnetic field environments of zero, 0.35 G 
transverse, and 0.35 G axial. The data associated with the zero and 0.35 G axial fields 
are very similar, but the 0.35 G transverse field causes the contactor plasma clouds to 
extend further downstream as shown in Fig. 18. 

The effects of magnetic field on the contactor plasma cloud and double layer 
regions at 1 G axial and transverse conditions can be seen by comparing the data of 
Figs. 19 and 20. Figure 19 shows the contactor plasma cloud extends further 
downstream and the double layer voltage drop increases from zero to 55 V as the 
electron collection current is increased from 50 to 200 mA in a 1 G axial field 
environment. This observation is in agreement with the unpublished results of previous 
electron collection experiments conducted at CSU and LeRC when low contactor 
flowrates were used and no ignited mode transition was observed. Results similar to 
those shown in Fig. 19 were also obtained when the magnetic field was set to zero. 

Figure 20 shows plasma potential profiles obtained at a 1 G transverse 
magnetic field condition. The profile corresponding to 50 mA of electron collection 
displays one double layer. However, as the electron collection current is increased to 
100 and 150 mA, two double layers develop. In addition to multiple double layers 
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Fig. 19 Effect of Electron Collection Current on Plasma Potential Profiles in a 1.0 G Axial 
Magnetic Field 





occurring at higher electron collection currents, higher plasma noise was observed. 
Unfortunately, the noise data were qualitative and numerical comparisons cannot be 
made. The noise level tended to be greatest through the multiple double layer regions 
and much lower close to the contactor and in the ambient plasma region. 

Up to this point, emphasis has been placed on the axial variation of plasma 
properties. The uniformity of the ambient plasma properties in the radial direction are 
shown in Fig. 21. This figure contains plasma potential, density and temperature data 
at the nulled, 0.35 G transverse, and 0.35 G axial magnetic field conditions at fixed 
axial positions of 61, 61, and 43 cm, respectively. The plasma potential is shown to be 
relatively constant at about 6 V over the 1 m radial region investigated, although it 
increases to about 8 V at a radial position of -50. It is noted that the negative radial 
positions correspond to locations between the centerline and the cryo-pump side of the 
facility (see Fig. 13) and positive radial positions correspond to locations between the 
centerline and the quadropole mass analyzer. Both negative and positive radial position 
data were shown in Fig. 21 in order to display the level of symmetry in the ambient 
plasma. The middle plot in Fig. 21 contains plasma density data, which show 
maximum densities on the centerline. Although the plasma density data show some 
scatter, the density is generally highest for the axial magnetic field condition and lowest 
for the transverse magnetic field. The electron temperature is nearly constant in the 
radial direction at about 1.2 eV as shown in the bottom plot of Fig. 21. 

The maximum magnetic field that could be induced in the IFSI facility was 1.6 
G, and plasma potential profiles corresponding to a transverse field of this magnitude 
are shown in Fig. 22. The low electron collection current of 50 mA shows a single 
well-defined double layer which develops between 35 and 45 cm. At higher currents of 
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Fig. 22 Effect of Electron Collection Current on Plasma Potential Profiles in a Strong 1.6 G 
Transverse Magnetic Field 




100, 150, and 200 mA, two double layers are shown to develop which are less well- 
defined and extend further downstream. In addition, as the current is increased, the 
total voltage drop across the double layer increases. 

Conclusions 

Similar double layers exhibiting large voltage drops and very distinct 
boundaries are observed in electron collection experiments conducted in both the IFSI 
and CSU laboratories. Independent plasma potential measurements made using 
Langmuir and emissive probes in regions away from the double layer agree well 
(typically within 1 V). However, emissive probes are attractive because they yield 
potential data that vary continuously through the double layer while Langmuir probes 
yield reliable plasma potential data only on either side of the double layer. Further, 
emissive probe results facilitate interpretation of the Langmuir probe data. A 
comparison of CSU and IFSI Langmuir probe analysis procedures suggest that both 
have advantages and disadvantages. The CSU analysis procedure can be used to model 
high energy, non-thermal electrons that are frequently present in the contactor plasma 
cloud. The IFSI procedure is easier to use when nearly-Maxwellian plasmas are being 
investigated and plasma potentials are being determined from thick-sheath Langmuir 
probe traces. 

The effects of magnetic field strength (0 to 1.6 G) are small when the field is 
oriented along the direction of current flow between the contactor and simulator. 
Increases in the strength of a magnetic field oriented transverse to the current flow 
direction induce multiple double-layers and increased noise levels near the double 
layers. 
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Appendix A 


SIMPLE TECHNIQUE FOR OBTAINING ELECTRON ENERGY DISTRIBUTIONS 

Introduction 

In order to estimate the electron energy distribution function directly in a 
relatively low density, isotropic plasma using a spherical Langmuir probe, it is 
necessary to compute the second derivative of its current/ voltage characteristic curve. 

It is generally very difficult to differentiate experimental data twice without amplifying 
the inherent noise in it to the point where it dominates any useful information. 
However, when clean, smooth experimental data are obtained using a device with a 
very low-pass filter and averaging capabilities like those available on Keithley 617 or 
237 programmable electrometers and special numerical procedures are performed, 
realistic electron energy distributions can be generated. This appendix presents 
information about a very simple numerical method, which can yield convergent, 
relatively smooth derivatives of experimental data. 

Application Example and Numerical Approach 
A Langmuir probe trace typically consists of many discrete current/ voltage data 
pairs [V n , J(V n ) - n = 1, 2, ... , N] (equally-spaced in voltage). Figure A1 shows 
two typical, thick-sheath Langmuir probe traces constructed from plotting discrete 
current/voltage data sets. The probe voltage for these traces is referenced to the tank 
ground of the CSU facility, and they were measured in the contactor plasma cloud and 
ambient plasmas during a test in which a 50 mA electron current was being collected 
by the contactor from the ambient plasma. The trace obtained in the contactor plasma 
cloud contains features which suggest that a low-energy group (probably Maxwellian) 
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Fig. A1 Typical Langmuir Probe Current/ Voltage Traces Measured in the Ambient and 
Contactor Plasma Cloud Regions and Corresponding Second Derivatives 



and a higher energy group of electrons are present, while the trace corresponding to the 
ambient plasma appears to contain only one, low-energy group. These characteristics 
can be seen more clearly by examining the second derivative curves shown beneath the 
current voltage traces. The second derivative curves can be used along with plasma 
potential measured by an emissive probe to compute the electron energy distribution 
function [18] using the equation 


F(E) 


? 3/2 1/2 j 

2 m ‘ jE dJ 


e 3 , 2 A, 


dV L 


. (Al) 


In Eq. (Al), E is defined as the difference between plasma potential and any given 
probe potential (V). 

In order to obtain the second derivative required in Eq. (Al), the discrete data 
points in a Langmuir probe trace were modified using the following procedure. First, 
a straight line, which connects the two end points of the trace, was subtracted from the 
data set: 


In Eq. A2, x and y represent potential and modified current, respectively. The x and y 
data pairs are next represented as a continuous function (a Fourier sine series which 
was found using a least-squares fit) 


i.e. 


y 


E A j sin 

;= i 


2rjx 


. (A3) 
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The number of terms in the series (m) was chosen to be half of the number of points in 
the data set (N) in order to avoid aliasing. Finally, the coefficients of the sine series 
were multiplied by the following convergence factor given by Lanczos [15] 




sin <^> 


LL 

N 


• (A4) 


The new sine series composed of the Bj coefficients can be differentiated analytically. 
Finally, the overall procedure can be repeated to obtain the second derivative required 
in Eq. (Al). 

Figure A2 contains two electron energy distribution functions which correspond to 
the Langmuir probe traces shown in Fig. Al that were normalized to their most 
probable value. When the Langmuir probe traces were analyzed using traditional 
techniques, they exhibited a temperature of about 3 to 4 eV for the low-energy 
(Maxwellian) electron groups. However, the most probable energy of both 
distributions is about 3.5 eV and this value is higher than expected (if the low-energy 
group electrons were Maxwellian, this result suggests that their temperature would be 
about 2 * 3.5 eV = 7 eV). Another measure of a distribution is its full-width, half- 
maximum (FWHM) value. The distributions shown in Fig. A2 display FWHM values 
of 7 to 8 eV which correspond to Maxwellian distribution temperatures of 3.9 to 4.5 
eV (i.e. T e ~ FWHM/1.8). This result agrees rather well with the Langmuir probe 
analysis estimates of 3 and 4 eV. 
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ELECTRON ENERGY [E] (eV) 

b. Ambient Plasma Data 


Fig. A2 Electron Energy Distributions Calculated from the Langmuir Probe Traces 
Shown in Fig. Ai 




Conclusions 

The distribution functions obtained by computing the second derivative of 


experimental Langmuir probe traces are typically too noisy to be of any value, 
however, they can be made less noisy through the use of a Lanczos convergence factor 
in conjunction with a sine series fit. The use of the Lanczos convergence factor is 
equivalent to smoothing the experimental data, and it causes smoothing errors. In 
addition, errors caused by inaccuracies in plasma potential measurements and natural 
rounding of the Langmuir probe trace near plasma potential in a noisy plasma probably 
introduce some inaccuracies into the electron distributions functions. Although errors 
reduce the accuracy of the computed distribution functions, the procedure outlined here 
does provide useful, qualitative estimates of electron energy distributions. 
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Appendix B 


PLASMA POTENTIAL MEASUREMENTS IN LOW 
DENSITY PLASMAS USING AN EMISSIVE PROBE 

Experiments were conducted at CSU, LeRC and IFSI-CNR to investigate the 
accuracy of plasma potentials measured in low density plasmas using emissive probes. 
The LeRC vacuum system in which the most extensive tests were conducted and the 
electrical circuits that were used to power the probe and plasma source system are 
shown schematically in Fig. Bl. Typical tests were initiated by striking a discharge 
within the ion source and allowing the plasma which it produces to expand into the 
chamber. The characteristics of the plasma were studied with an emissive probe that 
was fixed at radial and axial coordinates of —0.7 m and — 1.2 m, respectively, from 
axes defined by the ion source/facility centerline and the source aperture. 

The ion source is a simple, divergent field discharge chamber which is equipped 
with a hot filament cathode. One side of the cathode was connected to the tank wall 
and the ion source body (which also serves as the anode) was biased 50 V positive of 
the cathode using the anode supply shown in Fig. Bl. The discharge current (J A ) 
flowing from the cathode to the body was controlled by adjusting the cathode filament 
temperature using the heater supply (i.e. by adjusting J DF )- The source was operated 
on argon for all tests. 

The emissive probe is constructed from tungsten wire (76 /xm dia. and - 1.2 cm 
long) formed into a semi-loop with its ends attached to lead wires via low resistance 
contacts which are insulated from the plasma using ceramic adhesive. The high power 
potentiometer shown in Fig. Bl could be used to adjust the heating current (Jjj) over a 
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EMISSIVE PROBE CIRCUIT 


Fig. B1 NASA Lewis Electrical and Mechanical Schematics 








range from 0.5 A to 1.2 A. At the lower limit (0.5 A) the filament appeared dark to 
the naked eye, while it appeared white hot at currents over 0.9 A. A typical 
incandescent light bulb filament temperature is about 2800 K and this temperature is 
comparable to the emissive probe temperature at the higher heating currents. The 
potential of the common point between the precision 7 kfi resistors could be varied 
using a Keithley 237 electrometer and voltage source, and the electron current collected 
by the emissive probe could be measured as a function of this potential. Note that the 
Keithley 237 power supply sees an input impedance due to the two precision resistors 
and the effective impedance between the emissive probe and the plasma in which it is 
immersed. When the emissive probe is cold and is in a relatively low density, low 
temperature plasma, the effective impedance between the probe and the plasma is very 
high compared to 3.5 kQ— the effective impedance to the probe through the precision 
resistors. Consequently, the effect of the precision resistors on the current/voltage 
characteristic is small and can be neglected. However, when the emissive probe is 
heated white hot and is biased negative with respect to the plasma, its effective 
impedance to the plasma decreases. At this condition the precision resistors can affect 
the probe current/voltage characteristics by limiting the electron emission current to the 
surrounding plasma. 

Several current/voltage characteristic curves obtained using the emissive probe are 
shown in the upper plot of Fig. B2. They were measured at heating currents (J H ) of 
0.5, 0.8, 0.9, and 1.1 A when the ion source was being operated on argon at an anode 
current and voltage of 0.8 A and 50 V, respectively. When the probe was relatively 
cold (J H * 0.5 A), it behaved like a conventional Langmuir probe and the associated 
probe current/voltage curve could be analyzed to obtain a plasma potential of 
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3.8 V, a plasma density of 2.1xl0 6 cm* 3 , and an electron temperature of 3.1 eV. The 
pressure in the facility during this test (P 0 = 2.2xl0' 5 Torr) was determined by 
multiplying the value read from an ionization pressure gauge by the appropriate 
correction factor (0.66 for argon). 

As an emissive probe is heated to higher temperatures and held below plasma 
potential, it emits a progressively greater electron current. For example, the upper plot 
in Fig. B2 shows that increasing the heating current from 0.8. to 1.1 A causes the probe 
current measured below plasma potential to change from zero to large negative values. 
The bottom plot in Fig. B2 was constructed by subtracting the curve labelled J H = 1.1 
A from the one labelled 0.5 A. This plot shows that the emission current from the 
white hot probe increases as the probe is held more negative of plasma potential and 
that it does not saturate. Previous experiments conducted at CSU in higher density 
plasmas showed that the emission current does saturate and the reason it did not 
saturate in the Lewis tests is not understood. Possible explanations for this difference 
include a) higher neutral densities in the Lewis facility which enhanced ion production 
near the probe and facilitated increased ion collection by the probe as it was biased 
more negative and b) limitations imposed by the 7 kQ precision resistors on the 
emission capabilities of the hot probe. The effective impedance associated with 
operation of the hot probe at potentials below about 0 V in Fig. B2 is simply the 
reciprocal of the absolute value of the slope over this same potential region. This turns 
out to be —20 kfi, which is about an order of magnitude greater than the effective 
impedance associated with the precision resistors (3.5 kfl). This result indicates that 
the effective impedance of the precision resistors should have a negligible effect on the 
current\voltage characteristic curve. It can be concluded from the data of Fig. B2 that 
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the actual emission characteristics of a hot emissive probe are highly dependent on the 
probe temperature and the properties of the plasma in which it is immersed. A model 
that describes probe emission characteristics like those in Fig. B2 would have to 
account for both of these effects. 

The characteristic curves shown in the top plot of Fig. B2 also differ at potentials 
above plasma potential— namely, the hotter probes collect greater electron currents. It 
is noted that the emissive probe temperature at the highest heating current is only about 
3000 K so it should be unable to emit electrons when it is at any potential greater than 
about 0.3 volts above plasma potential. Consequently, all of the curves in the upper 
plot of Fig. B2 should fall on top of one another at potentials about 0.3 V above 
plasma potential. Several phenomena could induce the systematic probe temperature- 
related increase in electron current to the filament from the plasma at potentials beyond 
this value. Specifically, increases in probe temperature could induce: a) physical 
growth of the probe surface area, b) growth of the effective probe collection area as a 
result of increased conductivity of the insulating surface in contact with the probe wire 
and c) increased rejection of contaminants from the probe surface. Detailed evaluation 
of each of these possibilities suggests: 

a) The linear thermal expansion coefficient for tungsten when it is heated from 
293 K to 3000 K is about 1.6% [18]. If the tungsten wire expands in both length 
and diameter by this amount the surface area only increases by 2.7%. 

Conclusion- effect is negligible. 

b) The ceramic adhesive used at CSU is Ceramabond 569 and it has been found 
to become increasingly conductive as its temperature is increased above — 400 K. 
Depending on how quickly heat is conducted away from the ceramic adhesive in 
contact with the emissive probe, it is possible that some ceramic will become 
sufficiently conductive to increase the effective area of the probe exposed to the 
plasma as probe temperature is increased. Conclusion- effect is a possible 
contributor to observed error. 
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c) Contamination was observed on large Langmuir probes which were located in 
the facility during the three days of testing. However, analysis of the cold 
characteristic curve for the emissive probe indicated a plasma potential close to 
the point where the hotter curves begin to display electron emission. This 
agreement suggests that the contamination layer resistance at cold probe 
temperatures is not significant compared to the impedance between the probe and 
the plasma. Conclusion- effect is not considered to be substantial for the low 
density and temperature plasmas being investigated. 

It is also possible that plasma conditions varied as probe heating current was changed 

but such changes would be expected to be random rather than systematic. The data in 

the upper plot of Fig. B2 show a systematic change in probe current with heating 

current so this is considered an unlikely explanation. 

Figure B3 is a plot of floating potential versus emissive probe heating current 

corresponding to the data shown in Fig. B2. The floating potentials were measured by 

noting the potential where the current sensed by the Keithley 237 was zero and they are 

shown to remain at about -10.7 V until the heating current is about 0.7 A. At heating 

currents above 0.7 A, the floating potential rises quickly and appears to saturate at a 

value between 4 and 4.6 V. This number is in reasonable agreement with the value 

obtained for the plasma potential from the Langmuir probe analysis program (3.8 V). 

Consequently, Fig. B3 suggests that the emissive probe is performing properly and that 

its floating potential at high heating currents is a good indication of plasma potential. 

In order to change the plasma density at the emissive probe location, the ion 

source discharge (anode) power was varied. Figure B4 shows the effects of changes in 

anode current (J A ) on the plasma density and plasma potential sensed by the emissive 

probe. One important observation that can be made from the data shown at the top of 

Fig. B4 is that the plasma density varies linearly with anode current. Note that the 

anode voltage (V A ) was held constant at 50 V so discharge power should be directly 
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Fig. B3 Effect of Heating Current on Emissive Probe Floating Potential 
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Fig. B4 Application of Emissive Probe: Effect of Anode Current on Plasma Density and 
Potential 
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proportional to the anode current for these curves. This suggests that plasma density in 
the vicinity of the emissive probe is linearly proportional to ion source discharge 
power. The electron temperature was found to be relatively constant at about 3 eV 
over the range of anode currents shown in Fig. B4. The bottom plot in Fig. B4 shows 
that the plasma potential was also relatively constant at about 4 V . 

The plasma produced within the ion source should have a potential near the 
source anode potential and electrons and ions created in the source should be drawn 
from this potential into the facility. The plasma potential in the facility should adjust 
itself so that the rates of ion and electron loss from it will match the rates of ion and 
electron supply from the source. This potential should be negative of the source 
plasma and positive of the tank wall so electric fields in sheaths near the tank wall can 
allow most ions to flow to and most electrons to be repelled from the tank walls. 

Under these conditions ions would be accelerated and energetic electrons would be 
decelerated (and cooled) as they travel from the interior of the source to the facility. 
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